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ABSTRACT 

We present the global distribution of fine structure infrared line emission in the Cassiopeia A 
supernova remnant using data from the Spitzer Space Telescope's Infrared Spectrograph. We identify 
emission from ejecta materials in the interior, prior to their encounter with the reverse shock, as 
well as from the post-shock bright ring. The global electron density increases by > 100 at the shock 
to ~ lO^cm"^, providing evidence for strong radiative cooling. There is also a dramatic change 
in ionization state at the shock, with the fading of emission from low ionization interior species 
like [Sill] giving way to [Siv] and, at even further distances, high-energy X-rays from hydrogenic 
silicon. Two compact, crescent-shaped clumps with highly enhanced neon abundance are arranged 
symmetrically around the central neutron star. These neon crescents are very closely aligned with 
the "kick" direction of the compact object from the remnant's expansion center, tracing a new axis 
of explosion asymmetry. They indicate that much of the apparent macroscopic elemental mixing 
may arise from different compositional layers of ejecta now passing through the reverse shock along 
different directions. 

Subject headings: supernova remnants - infrared: general - x-ray: ISM - supernovae: individual 
(Cassiopeia A) 



1. INTRODUCTION 

Supernova remnant Cassiopei a A is both n earby 
(3.4kpc) and very young (^330yr,"Fese n et al.|2006 ), giv- 
ing it a bright, richly detailed ejecta structure which has 
led to intensiv e study at many wavelengths, from gamma 



rays to radio CRudnick| |2QQ2| [Albert et al.||2007D . Spec- 
troscopy of distant light echoes from the origmal blast 
indicates that Gas A was the type lib core-collapse of a 
~15M0 main sequence star (Krause et al. 2008), which 
models indicate are highly stratified, having lost a signif- 
icant portion of the original outer layers o f hydrogen and 
helium prior to collapse in stellar winds (Woosley et al. 
1987). An initial outgoing "forward" shock wave driven 
by th e blast is seen as a thin X-ray edge (Gotthelf et al. 
2001) expanding at ~5000km/s, with a reverse shock 
being driven back into the outgoing ejecta as it interacts 
with circumstellar material and the interstellar medium, 
expan ding at roughly half the rate of the forward shock 
(jDeLaney & Ru dnick|[2QQ3 ). 

The high resolution available (l^'=0.016pc) makes Gas 
A ideal for tracking the physical conditions in the lay- 
ered ejecta material, following nucleosynthesis both be- 
fore and during the supernova explosion, the subsequent 
mixing and inhomogenous transport of the ejected mate- 
rials, and the formation of dust in situ. The mid-infrared 
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(MIR) is an ideal wavelength for studying these pro 
cesses, as it is minimally affected by the significant line of 
sight dust extinction {Ay = 4.6 - 6.2, Hurford & Fesen 



1996 ), and probes a rich suite of abundant ionized species 
including Si, S, Ne, and Ar which range in ionization 
energies from 8 - lOOeV. Previous MIR studies of the 
remnant using Infrared Space Observatory (ISO) map- 
ping placed constraints on the dust formation and the 
appearance of macroscopic mix ing of materials (|Lagage 
et al.|1996 Douvion et al.|19'99 ), and explored the veloc 



ity shifted line and dust emission in various regions using 
high-resolution MIR spectroscopy ( Arendt et al.||1999p . 

We report on the first large, sensitive infrared spec- 
tral maps of the remnant, obta ined by the Spitzer Space 
Telescope (Werner et al. 2004). The distribution of in- 
terstellar material, ejecta, and shocked circumstella r ma - 
terial is shown in the 24 /im maps of [Hines et al.| ( |2QQ4[ ). 
A model of the composition and mass of emitting dust 
associated wit h the remnant us i ng this MIR data set was 



undertaken by Rho et al.|( |2QQ8| ). |Ennis et al.|(|2QQ6j ) high- 
lighted the associated deep 3.6-8 um Spitzer IR AC imas;- 



lighted the associated deep 3.6-8 //m Spitzer IRAC imag 
ing, exploring the survival of the various nuclear burning 
layers. Here we present the distribution of the global fine 
structure infrared emission lines and constrain conditions 
in the emitting environments over the full remnant. The 
complex Doppler velocity field as trace d by these fine 
structure lines will be detailed separately (DeLaney et al. 
.2008^ . 

2. OBSERVATIONS AND REDUCTION 

Low-resolution 5-38 /im spectral maps were obtained 
in January, 2005 w ith the Infrared Spectograph (IRS, 
Houck et al. 2004) aboard Spitzer. A single mapping 



observation consisting of 364 individual rastered spectra 
was obtained in the Long-Low module (LL, 15-38 /im), 
utilizing 5.08'^ steps across the slit (1/2 slit width), and 
158'^ steps along the slit (approximately the full sub-slit 
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length), for a total of 12.6 seconds per position. Short- 
Low (SL, 5-15 /im) coverage was obtained by tiling four 
separate spectral maps to accomodate the maximum ob- 
servation duration limits of the IRS, comprising 1218 
total independent spectra. Precise scheduling was re- 
quired to fully tile the remnant without allowing gaps to 
open between the quadrants. The SL slit was stepped 
at 2''x52'', for slightly less than double coverage and an 
effective exposure time of 11.4s per position. The best 
effective resolution is ^2" in SL and ~8'' in LL (varying 
linearly with wavelength). The total angular coverages 
obtained were llfOxTfS (LL) and 6:3x5.^9 (SL). 

Background removal was performed in situ^ using spec- 
tral pointings selected from the periphery of the maps, 
choosing regions with low relative 24 /im surface bright- 
ness. Since the IRS sub-slits (mapped to separate grating 
orders) are separated by 180^' (LL) and 68^' (SL), multi- 
ple positions outside the main remant are obtained in the 
course of a full spectral map in both orders, which aids 
in selecting suitable background spectra. However, since 
Cas A lies close to the Galactic plane, foreground and 
background emission from intervening and background 
Galactic cirrus clouds — including molecular hydrogen, 
PAH emission, and weak low ionization emission lines — 
is present throughout the region. This complex filamen- 
tary emissi on, readily seen in extended 24 /am maps ( Kim 
et al.|2QQ8| ), varies considerably over the face of the rem- 
nant, and is also present in the selected background re- 
gions. Although the contaminating line emission is very 
faint — ^2% of the line surface brightness of the bright 
emitting ring — it provides a modest limitation on our 
ability to recover fluxes for the lowest surface brightness 
lines (see also Fig. [T]). 

All IRS data were reduced from pipeline version S17 
data products with Cubism, a custom tool created for 
the assembly and analysis of spatially-resolved spectral 
cubes from IRS spectral maps ( [Smith et al.| |2QQ7). A 
total of four separate spectral cubes were created, one 
each for orders 1 and 2 of both SL and LL. The flux cal- 
ibration was achieved using appropriate extended-source 
FLUXCON corrections, and aberrant or "rogue" pixels 
were automatically cleaned using the statistical methods 
built into Cubism, with a small amount of additional 
manual removal. The short 6.3s exposure ramps contain 
only four samples, which leads to occasional anamolous 
cosmic ray correction, inducing small bad pixel clusters. 
These were mitigated with statistical record-level bad 
pixel flagging. 

3. SPECTRAL MAPS 

A full 7.5-14.7 /im map created from SL order 1 is 
shown in Fig. [l] and o ffers excellen t correspondence with 
the 8 /im IRAC map of Ennis et al. ( [2006} . Two represen- 
tative full spectra extracted from the interior and along 
the bright northern rim are also shown. These two spec- 
tra highlight the large dynamic range of line intensity, 
line ratios, and continuum shape across the remnant face. 
The broad peaks at 9 and 21 /im aris e from warm pro to- 
silicate and other dust emission (see |Rho et al.|[2QQ8l ) . 

Surface brightness maps from the eleven brightest 
emission lines (readily seen in Fig. [T]) are shown in 
Fig. [2] Spectral line maps were created by integrat- 
ing the spectral cube intensities over a wavelength re- 
gion large enough to sample the line flux at all veloci- 



TABLE 1 
Detected Emission Lines 



Line 
(1) 


Wavelength 
(2) 


Peak Intensity 
lO-'^W/m^/sr 
(3) 


Total Flux 
10-13 w/m2 
(4) 


Niii] 


6.636 


2.4 ±0.85 




At ii] 


6.985 


221. ±0.86 


7.72 


Arm] 


8.991 


41.7±0.82 


1.24 


IV] 


10.511 


21.0 ±0.43 


1.17 


Ne 11] 


12.814 


33.7±0.21 


1.41 


;Nev] 


14.322 


1.75±0.17 




^Neiii] 


15.555 


8.62 ±0.17 


0.49 


Fe 11] 


17.936 


2.49 ±0.14 




III] 


18.713 


13.2 ±0.23 


1.44 


^0 iv]^ 


25.890 


27.1 ±0.29 


4.95 


Peii]^ 


25.988 






III] 


33.481 


4.08 ±0.18 


1.02 


Sin] 


34.815 


12.3 ±0.21 


3.16 



Note. — Col. (1): line name. Col. (2): rest wavelength. 
Col. (3): peak surface brightness over 1.85''^ (gL) g^^d 5.08''^ 
(LL) areas. Col (4): total flux integrated over the remnant in a 
5.^3x5.^3 region; statistical uncertainties in the total are 0.01 or 
less in these units. Faint lines excluded. 

^ [O iv] and [Fe ll] are blended into a single emission peak at LL 
resolution. 



ties present within the remnant (roughly -5000 km/s to 
+7000 km/s, see DeLan ey et al. 2008) . Prior to integra- 
tion, the local dust and synchrotron continuum emission 
was removed by subtracting a weighted average of flank- 
ing wavelength regions, carefully chosen to avoid con- 
tamination by nearby lines at all velocities. The aver- 
aging weight used at each wavelength within the line re- 
gion {Xiine) for a given continuum wavelength (Xcont) was 

I'^cont Xiiyiq I 

The broad dust feature present throughout most of 
the ejecta has a sharp peak at 21 /im, which induces 
a steep and variable continuum in the 15-21 jiia range 
(see Fig. [T]). Estimating accurate line map fluxes in 
this range (LL order 2) required a preliminary contin- 
uum removal step. The cube was flrst smoothed spa- 
tially by three pixels. A linear or quadratic background 
function at each (smoothed) spatial position was fltted 
using two to four wavelength intervals free of line con- 
tamination, separately chosen for each line in the range. 
This noiseless background estimate was removed from 
the unsmoothed cube, and flanking continuum regions 
were then subtracted as above to eliminate any resid- 
ual continuum contributions underlying the line. For the 
faint [Fell] 17.94 /im line in this region, uncertainty in 
the continuum subtraction dominates the residual sys- 
tematic errors. 

Table [l] lists the lines detected at low resolution, with 
the peak surface brightness along the bright ring, as well 
as integrated line flux over the entire remnant. Very faint 
lines known to be associated primarily with Galactic cir- 
rus emission, including rotational H2 S(0)-S(2) and PAH 
emission features, are omitted, as are the integrated flux 
measurement of the faintest ejecta lines, which would 
be dominated by systematic errors. Among the brighter 
lines, systematic uncertainties in the total flux are 10% 
or less. Signal to noise in the brightest regions of the 
line maps was typically 30 or above. [Am] is by far 
the brightest line detected, and as it also offers the best 
spatial resolution, has proven very useful for Doppler de- 
composition of the remnant. 
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Fig. 1. — At left, a full 7.5 /im-14.7/im map created from the SLl cube, 320^^x320^^, square root scaled, with offset positions labelled 
relative to the expansion center (23:23:27.77, 58:48:49.4; Thorstensen et al. 2001). Representative spectra extracted in circular apertures 
are shown at right. Lines detected from the ejecta are labeled on the spectrum ot the bright northern rim at top, and residual PAH features 
associated with the Galactic foreground are indicated on the much fainter interior emission spectrum below. Broad peaks at 9 and 21 fim 
are due to dust emission (see Rho e t al.,2008) , and the periodic signal at 20-25 /xm is instrumental fringing. The wavelength range used to 
construct the map is indicated with vertical dotted lines. 



In many cases line identification is hampered by 
Doppler broadening and, less frequently, multiple line-of- 
sight velocity components. The [O iv] line at 25.89 jam 
is a special example of such a case, as it is profoundly 
blended with [Fell] 25.99 /im even at high spectral res- 
olution. The line map at 26 /im therefore contains con- 
tributions from both of these lines. Constraints on the 
fractional contribution of [Fell] and [Oiv] in different 
regions can be made using expected line ratios to other 
lines of Fe, inluding [Fell] 17.94 /im, as well as [Felll] 
22.93 /im, detected in our follow-on high-resolution maps 
of selected regions. Such an analysis will be considered 
in a subsequent study, but is consistent with both [Fe ll] 
and [Oiv] contributing and in differing amounts across 
the remnant. 

In cases where Doppler confusion could support more 
than one line identification (most significant at the rem- 
nant center), we used our high-resolution spectra to re- 
solve the ambiguity. An example for [Sill] is given in 
Fig. [3) A consistent set of three velocity components is 
visible for [S ill] , [Si ll] , and a faint [Fe ll] line undetected 
in our low-resolution data set. 

A bright ring of emission, with overlapping smaller sub- 
rings most clearly seen in the N and SW, is the dominant 
structure evident in all line maps, but significant varia- 
tions in the line strengths occur throughout, in particular 
when comparing the interior to the ring. For the most 
part, the line strength variations on the bright ring are 
not complicated by the three-dimensional structure of 
Cas A. This is apparen t in the Doppler reconstruction of 
Lawrence et al. (1995) which shows that there are very 
few overlapping structures along the same line of sight, 
with the exception of the NNE and the base of the NE 
jet. The same cannot be said of the interior emission, 
where nearly every line of sight has two or mor e domi- 
nant Doppler components ( DeLaney et al.|[2QQ8 ). 



4. PROGRESSION OF THE SHOCK 

A reverse shock driven by the interaction of high ve- 
locity ejecta with the circumstellar material is propagat- 
ing backwards into the outflowing material at a speed 
of 250 0-4400 km/ s (relativ e to the bulk outflo w velocity, 
DeLa ney fc Rudnick|[2QQ3| [Morse et al.||2004D . Material 
inside the reverse shock boundary was "once- shocked" 
by the initial outgoing blast wave, but has subsequently 
cooled and recombined, and has not yet been re-heated 
and re-ionized by the returning reverse shock. Three 
dimensional Doppler reconstructions of optical and in- 
frared lines such as [Am] show that they illuminate the 
reverse shock in a broad band appr oximately in the plane 
of the sky, forming t he bright ring ( Lawrence et al.|1995 



DeLaney et al.|2QQ8 ). There is little or no emission from. 



these ring lines across the face of the remnant; how- 
ever, interior to the ring, strong emission lines of the 
low ionization species of S and Si appear. The veloc- 
ity distribution of the interior [Sill] emission, shown in 
Fig. [4j demonstrates that this central material lies fully 
inside the reverse shock boundary. In addition to these 
lines, strong 26 /im emission is present in the interior, 
contributed by [Fell] and/or [Oiv]. 

4.1. Electron Density 

We can probe the change in density of the ejecta as 
they encounter the reverse shock using the [S ill] lines. A 
temperature-insensitive density diagnostic is available in 
the 18.71 /im/33. 48 /im ground-state fine structure dou- 
blet line ratio, arising from low- lying collisionally excited 
levels. Recovering electron density requires balancing 
the rates of collisional excitation and de-excitation and 
radiative transitions into and out of the relevant levels 
of the ion. For ions with p^, p^ and p^ ground state 
electron configurations, like S++, it is sufficient to con- 
sider only the 5 lowest lying energy levels ( Oster brock fc] 
Ferland|[2006 ) . Up to date temperature-dependent colli- 
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Fig. 2. — Surface brightness maps in the eleven brightest emission hnes, sorted by order of increasing wavelength. Each image is 
300^^x300^^ with offset positions as in Fig. ^ The integrated, continuum-subtracted line surface brightness is square root scaled between 
low and high thresholds specified in brackets in each pane (in units of 10~^ W/m^/sr). 
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Fig. 3. — A portion of the high-resolution IRS spectrum of the 
central remnant, used to identify [Sill] emission at 35.3 /im. The 
[Slll] 33.48 /im line was used to fit three velocity components at 
-1549 km/s, 1175 km/s, and 4303 km/s. These three components 
are marked for [S ill] (solid lines), [Sill] (dotted), and a weak line 
of [Fell] not detected in our low resolution data set (dashed). 
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Fig. 4. — Velocity distribution of [Sin] emission interior to the 
bright ring. Hatched regions represent the expected velocities (- 
4520 km/s, +6061 km/s) and the velocity width (1500 km/s) of the 
reverse shock at the remnant's center ( iReed et al.|1995j >. 
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Fig. 5. — The [S lll] 18 /im/[S lll] 33 /im density diagnostic for 
three representative temperatures. 
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Fig. 6. — The electron density in the remnant, computed from the 
[S III] 18 /im/[S III] 33 /im ratio. Regions without sufficient signal in 
the ratio map are shown gray. The darkest regions, which include 
most of the interior emission, are in the low d ensity limit of Ue 
100 cm ~^. The expansion center position of [Thorstensen et al.| 
( |2001| > is marked as a small circle, with offsets relative to it. 



sional rates were taken from the IRON project ("Hummer 
et al.p 993) using its TIPBASE database interface. Ra- 
diative rates were obtained either from the IRON project, 
or, where available and reliable, from the NIST Atomic 
Spectra Database^. The five coupled linear rate equa- 
tions were solved numerically for all level populations, 
imposing the constraint that the population fractions 
sum to unity. This solution is performed for a given 
temperature and density, interpolating the temperature- 
dependent collisional rates over log(T). The resulting 
line ratio dependence on density is illustrated in Fig. [5] 
The [Slll] 18 /im/[S III] 33/im line ratio saturates in the 
high density limit at 11.9, when collisions are rapid 
enough to populate the levels in equilibrium according 
to their statistical weights, and at a low density limit of 
0.56, when radiative effects dominate. 

To recover electron density in the remnant, the two 
[S III] line maps were spatially registered and convolved 
to matching resolution using a custom convolution kernel 
constructe d from Spitzer STinyT im PSF models in the 
manner of Gordon et al. (2008). The matched images 
were then ratioed, and corrected for Ay =5 magnitudes 
of extinction u sing the interstellar extinction curve of 
Chiar & Tielens (2006) (which results in a modest 10% 
increase). All data in the ratio map with signal-to-noise 
below 2.25 were discarded, and map regions within la 
of the low or high-density ratio limits were set to the 
limiting densities. 

The [S III] 18 /im/[S ill] 33 /im ratio ranges from 0.5 near 
the center to a maximum of 7.5 on the bright rim. The 
resulting projected density map is shown in Fig. |6] This 
density is averaged over spatial scales of 10'^ (0.16pc). 
Although the bright ring is not generally confused by 
overl apping Doppler compo nents at this spatial resolu- 
tion (Lawrence et al. 19951), high-resolution optical and 
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X-ray images resolv e structure out at 0.2-1'^ (Fesen et al. 
2QQl| iHughes et al.||2QQQD , so even without projection ef- 



fects, some blending of structure is unavoidable. In re- 
gions where two or more strong emission components are 
present along the line of sight, which includes most of the 
interior emission, the [S ill] line ratio represents the av- 
erage emission- weighted electron density. However, in- 
dividual line ratios calculated for red- and blue-shifted 
[S III] line sets in the interior emission near the remnant 
center yielded consistent values of ~0.5. 

The interior [S ill] emission is uniformly in or near 
the low limiting density of the [Slll] diagnostic (ng ^ 
100 cm~^), whereas the electron density in gas com- 
pressed by the reverse shock is much higher, peaking 
above 1.5xl0^cm~^. This general result holds even as- 
suming much lower temperatures, down to 1000 K. Un- 
shocked interior emission at much lower densities has 
also been ob served in the type la supernova remnant 
SN1006 (e.g. i Hamilton fc Fesen] [1988 ), and were origi- 
nally identified in Cas A through free-free absorp tion of 
low frequency radio em ission by Kassim et al. ( p!995| ). 
Hurford fc Fesen (1996) obtain densities ranging from 
5000-15000 cm~"^ for a range of optical knot com posi- 
tions, consistent with |Chevalier & Kirshner (1979), and 
in good agreement with our results on the bright rim. 

Since the density enhancement in strong non-radiative 
shocks is limited to a factor of four, the strong density en- 
hancement provides evidence for rapid radiative cooling 
of the infrared emitting clumps after they pass through 
the shock. Alternatively, hidden dense clumps may be 
entrained in the outfiowing material, which are illumi- 
nated only as they pass through the shock. This could 
then represent an extension to higher densities of the 
dumpi ness apparent in the more diffuse X-ray-emitting 
ejecta QMorse et ar][2Q04D . 



4.2. Ionization State 

The progression of the reverse shock into the ejecta 
can be observed not only in the electron density, which 
increases by at least two orders of magnitude, but also 
in the ionization state of the outfiowing atomic material 
as it encounters the reverse shock. 

Figure [7| illustrates this progression by comparing line 
maps from two infrared lines arising from ions of dif- 
fering ionization potential — [Sill] (8.2eV), and [Siv] 
(34.8eV) — together with a summed pair of high-energy 
X-ray helium-like and hydrogen-like K-alp ha resonance 
lines: Sixill (0.5keV), and Sixiv (2.4keV) (jHwangeTaL 
2000[ ). There is clear layering of ionization state, from 
low energy species in the interior, higher energies on the 
bright ring where optical emission is also seen, and very 
high energies traced by X-ray line emission extending be- 
yond the infrared-bright rim itself. Though line of sight 
projection of multiple structures complicates the view in 
the northern rim, it is evident that the observed layering 
is highly non-uniform as a function of azimuth as well as 
along the line of sight. 

In various positions, more or less of the 0-burning 
layer has encountered the reverse shock. In the direc- 
tion along our line of sight, very little of this material 
has reached the shock. These regions are dominated by 
low ionization lines, [Sill] in particular (see Fig.!?]), which 
lie at lower velocities than the shocked ring ofemission 
(Fig. [4]). In these directions, neither element is seen in 




-150 



50 

Aa (arcsec) 



150 



Fig. 7. — A color coded map in [Sin] 34.8 /im (red, ionization en- 
ergy 8.2eV), [Siv] 10.5 /im (green, ionization energy 34.8eV), and 
the summed X-ray emission of Si xill+Si xiv (hydrogen and hehum- 
Hke sihcon, blue, ionization energy 0.5-2.4keV). Square-root scal- 
ing is used for all three channels, and the coordinate scale follows 
Fig. [2] with a 320^^ image size. Immediately evident is the inside- 
out layering of ionization state, with low energy species in the un- 
shocked interior region, and higher energy species layered within 
the non-uniform shock region itself. 



higher ionization states in the infrared or X-ray. This 
low-ionization pattern also appears at different azimuths 
along the bright ring, indicative of positions where O- 
burning ejecta moving approximately in the plane of the 
sky have yet to encountered the reverse shock. This emis- 
sion is consistent with material inside the reverse shock 
which has been ionized by ambient Galactic stellar radi- 
ation and/or illumination from the rim of shocked ejecta. 

In other directions, moderate ionization state mate- 
rial traced by [Siv] is present, with weak or absent X- 
ray emission. Finally, there are regions where only the 
high ionization Si X-ray lines are seen, predominantly 
at large radii. These variations likely represent different 
stages in the post-shock histories of the inhomogeneous 
ejecta. The low density X-ray filaments and diffuse ma- 
terial are gradually heated and ionized as they penetrate 
the reverse shock, with a time history that depends on 
the s hape of the ejecta density profile (Laming & H wang] 
|2003[ ) . Higher density clumps of material can suffer a dif- 
ferent fate, reaching temperatures of only ~10^ K as they 
undergo rapid radiative cooling immediately post-shock, 
and are seen in [S iv] . 

The detailed histories of shock-heated material requires 
considerably more modeling, which incorporates varia- 
tions in ejecta density and its profile, ejecta velocities, 
and irregularities in the location of the reverse shock. 
It is not clear to what degree the density and ioniza- 
tion state of the material we now see passing through 
the shock front is determined by initial conditions of the 
ejecta — diffuse vs. clumpy, lower vs. higher velocity 
— and to what degree the heating, photo- and collisonal 
ionization, and subsequent cooling by line radiation re- 
distributes material between these phases. What is clear, 
however, is that the ejecta undergo a striking and rapid 
change in density and ionization level as they encounter 
the reverse shock. 
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Fig. 8. — Composite image of [Ann] (blue), [S iv] (green), and 
[Neil] (red), showing two pronounced neon-rich crescent-shaped 
regions to the N and S (square root scahng). Axes as Fig. [t] Also 
shown are the kinematic center of the remnant (magenta circle, 
[Xhorstensen et al. 2001) and X-ray localization of the remnant's 
compact object (rea square, 1" to the south). Green lines indicate 
the Icr range of the "kick vector" direction of t he compact o bject 
from the ejecta's expansion center: 169° ±8.4° ( |Fesen et al.|2 006). 
The two regions of enhanced neon abundance lie very close to this 
projected direction. Several much smaller neon-enhanced regions 
lie in the West along the X-ray jet direction. 

5. APPARENT INHOMOGENOUS MIXING AND 
THE NEON CRESCENTS 

The Gas A ejecta are by no means homogenous in 
composition or uniform in their layering. Using optical 
spectroscopy, 
abundance in 
ucts S, Ar 



Ghevalier fc Kirshner ( 1979 ) demonstrated 
lomogeneities in the oxygen burning prod- 
and Ga among the outlying fast moving 
knots, and suggested differing completeness levels of oxy- 
gen burning could explain these variations. Macroscopic 
mi xing of similar nucleo synthetic products was suggested 
by Douvion et al. (1999) using ISO MIR line ma ps of Ar, 
Ne, and S c overing a portion of the bright rim. [Hughes 
et al. ( |2QQQ ) used Ghandra X-ray spectroscopy to identify 
spatial inversion of iron and oxygen burning products. It 
would be fair to say that non-uniformity of composition 
is the rule rather than the exception in Gas A. 

In a stratified model in which different layers of ejecta 
material are coming into contact with the reverse shock 
front in different directions, the appearance of individ- 
ual emission lines depends on 1) the expansion velocity 
of material, 2) the position of the reverse shock with re- 
spect to it, and 3) the subsequent ionization state of the 
material as it passes through the shock. In this context, 
even if the gross nucleosynthetic layering is preserved in 
the outflow, with lighter elements overlying heavier ones, 
differing outflow velocities in different directions can gen- 
erate an apparent mixing of layers, since different strata 
are now arriving at and passing through the reverse shock 
boundary in different regions. Combined with line of 
sight projection effects, this can profoundly complicate 
interpretation of apparent elemental mixing, a conclu- 
sion also drawn in the jet-induced explosion model of 
Wheeler et aT] ( |2QQ8p . 

A powerful example of this apparent mixing is found 
in the neon maps. Neon is among the most abundant 



elements in the carbon and oxygen shell-bur ning cores 
of intermedi ate mass supernovae progenitors (Meakin & 
Arnett|2QQ6[ ) . In contrast to the optical neon lines, which 
are highly obscured by line of sight extinction, and there- 



fore difficul t to study (e.g. Ghevalier & Kirshner 1979 



Fesen| 1990), the MIR neon lines are bright and readily 
detected throughout the ejecta ring. 

A three color composite image of [Arm], [Siv], and 
[Nell] in Fig. [s] reveals two bright, crescent-shaped 
clumps of highly enhanced neon abundance to the North 
and South along the ring. Lines of Ar, Si, and S are 
very faint i n these regions, and they have excess 4.5 /im 
brightness ( [Ennis et al. 2006). The blended complex 
[O lv] + [Fe II] is strong as well, though the coarser angular 
resolution prevents completely separating the small re- 
gions from the surrounding ring material at these longer 
wavelength. 

The [Nell] /[Arm] ratio changes dramatically — from 
0.8 ± 0.4 on the rim outside of these crescent regions, to 
over 50 on the crescents themselves. Since the remnant 
regions occupied by the Ne crescents are not overly con- 
fused by overlapping Doppler structures (Lawrence et al. 
1995), the large ratio differences represent real abun- 
dance variations. Variations in abundance from optical 
emission lines of S and O a re also observed at spatial res- 
olutions of V ( Reed et'aL][l99 5 ) and down to i ndividual 



ejecta knots at^the 0V2 level ^besen et al.||2001 ). 

Ov erlaid on the image are th e remnant expansion cen- 
ter of Thorstensen et al. (2001) as well as the location of 
the X-ray compact object, which is offset from the kine- 
matic center by ^7'^ to the South. The projected direc- 
tion of motion, representing the asymme tric kick received 
by the compact object, is 169° ±8.4° ( ^Fesen et al.||2006D . 
The symmetric neon crescents lie nearly perfectly along 
this kick direction. Though the true-space velocity of rel- 
ative motion between the compact object and kinematic 
center is unknown, this strongly suggests a preferred axis 
of symmetry in the outflow along this direction. 

Of obvious interest is the total emitting mass of neon in 
the crescents. The observed line flux F of a collisionally 
excited emission line is: 



And? 



furiiondV 



(1) 



where Aui is the radiative rate between the two levels 
involved in the emission, fu is the fractional population 
of the upper level, nion is the ion density, V is the phys- 
ical volume over which the emission occurs, and d is the 
distance. Assuming the fractional density of the ion is 
constant over the emission volume, the emitting mass of 
the ion is therefore 



Mi, 



m 



lon't'ion 



Auifuhiy 



(2) 



where mion is the atomic mass. We take a representa- 
tive electron density on the bright rim of ng =4500 cm~^ 
from the [S III] analysis of § |4.1[ and bracket our uncer- 
tain ty in the temperature b y adopting T =5500-8000 K 
(see Hurford & Fesen 1996). Isolating the regions of ele- 
vated [Nell] /[Am] usmg polygonal apertures, we obtain 
the integrated fluxes for the crescents indicated in Ta- 
ble [| Each crescent is approximately 0.25 pc^ in size, 
and contributes roughly 10% of the remant's total [Nell] 
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line flux. The [Ne ill] line is also relatively strong in each 
crescent region, but [Nev] is very faint and not seen in 
the Southern crescent. The [Nev] /[Nell] ratio in the 
Northern crescent is > 8x that of the Southern region, 
although the velocity difference between the two lines 
at that position is > 1000 km/s, suggesting that the in- 
creased line of sight confusion in the Northern cap may 
play a role in this discrepancy. Ne^^ contributes negligi- 
bly to the mass. 

Assuming they emit in similar volumes, the densities of 
Ne^^ and Ne^^ relative to Ne^ can readily be estimated 
by forming a ratio of Eq. [l] for each of the higher ionized 
species to the lower, which eliminates all quantities ex- 
cept the line flux, upper level population, radiative rate, 
and frequency. We solve for the level populations using 
atomic data f or ne on obtained and applied in the manner 
described in §|4.1l This results in 77,^^++/?!^^+ = 0.17 in 



both regions, and, in the North, n^^A^/n^^^ = 0.003. 
Assuming these ions are the dominant source of elec- 
trons, and taking the fractional density of the unseen ion 
Ne^^^ (which does not emit in the MIR) as 0.05, we can 
estimate 



= 1 



nNe+++ 



. nNe4+ 



1.5 (3) 

Using Eq. |2] to calculate the mass, we find roughly 
1.8x10"^ Mq of neon in the crescent regions (dominated 
by Ne^ and Ne^^, and excluding neutral neon) . This 
mass is comparable to that found by Arendt et al. ( |1999[ ) 
in other (larger) segments of the rim. The entire ion- 
ized neon mass of the remnant, assuming an approxi- 
mately constant density throughout the bright rim, is 
8.6x10"^ M0, such that roughly 20% of the neon mass 
is concentrated in two small regions covering ^2% of the 
rim's area. 

The relative abundance of Ne to Ar (neglecting neu- 
tral species of both) can be calculated using line maps 
of [Am] and [Arm]. We assume the same temper- 
ature and densities used for neon, and calculate level 
populations as above. The relative abundance increases 
from nNe/^Ar=2.7 averaged over the full remnant to 
^Ne/^Ar=26 ou the S crescent region. These likely rep- 
resent lower limits on the full relative abundance, includ- 
ing neutral species, as the first ionization potential of Ne 
(21.6 eV) is substantially larger than that of Ar (15.8 eV), 
so that the average neutral fraction can be assumed to 
be larger in the former. 

Flowing outwards ~20 degrees from the plane of the 
sky at roug hly -5500 km/s (South ) and +4200 km/s 
(North, see DeLaney et al.||2008|), the ionized neon 
moving in these crescent-shaped clumps represent only 
4 X 10~^ of the esti mated 10^^ J of expansi on energy in 
the entire remnant (Willingale et al. 2003). Assuming 
the neutron star is moving in the plane 01 the sky, the 
7V0 separation it has obtained from the expansion center 
since the explosion (Fesen et al. 2006|) implies a veloc- 
ity of 342 km /s. Further assum ing it to have a mass of 
~1.4M0 (e.g.| Young et al.|2006 ), the momentum impulse 
imparted by the kick to the neutron star is then > 1000 x 
larger than that of either neon crescent. 

Since they represent such small fractions of the full 
eject a's kinetic energy and the momentum of the neu- 
tron star, it is unlikely that even a very large asymmetry 



in the neon crescents' initial velocities could have directly 
contributed to the kick experienced by the compact ob- 
ject with respect to the expanding ejecta. In many ways, 
however, the infrared emitting mass of neon in the cres- 
cents must be a significant underestimate of the full col- 
umn of material moving along this symmetric pair of 
directions, with layers of heavier elements which have 
not yet reached the reverse shock still largely unillumi- 
nated, neutral species missing from the estimate, hotter 
X-ray gas at greater radii, and other elements such as 
oxygen entrained in the crescents themselves not readily 
accounted for. 

In the final shell burning nucleosynthetic structure, the 
neon layer lies above the heavier si licon, sulfur, and ar- 
gon layers fWoosley & Weaver" 1995) . If the gross layering 
of pre-explosion elements is preserved, the strong over- 
abundance of neon in these crescents along the direction 
of the kick vector could in fact represent material which 
is late in meeting the reverse shock, having been accel- 
erated to slightly slower initial speeds. If so, it could be 
expected that the neon will gradually become less domi- 
nant as it progresses further beyond the shock and cools, 
with sulfur and argon filling in behind it. In this sense, 
this pair of directions is intermediate between the bright 
rim, cooling dominantly in lines of Ar, and the front and 
back surfaces along the line of sight, along which very 
little material has yet encountered the reverse shock. 

The fact that the compact neon crescents are nearly 
equal in size and in mass, and are arranged symmetrically 
about the expansion center, strongly indicates they arise 
not from macroscopic or turbulent mixing of the ejecta 
layers during expansion, but from a preferred axis in the 
initial outflow configuration along this direction. 

6. CONCLUSIONS 

We present global low-resolution 5-37 /im spectral 
maps of the Cassiopeia A supernova remnant obtained 
with Spitzer's IRS spectrograph. We detect numerous 
fine-structure emission lines arising from ions with ioniza- 
tion energies of 8-100 eV, and strengths that vary consid- 
erably across the remnant. The observed ions represent 
some of the most abundant species in oxygen and carbon 
shell burning supernova progenitor cores. We find: 

1. Strong interior emission of predominantly low ion- 
ization, heavy ionic species, dominated by [Sill]. 
These regions lie inside the bright rim along the line 
of sight, and probe slower moving material which 
has cooled from the passage of the initial outgoing 
forward shock wave, but which has not yet encoun- 
tered the reverse shock. 

2. The electron density in the remnant increases 
rapidly from < 100 cm~^ in the unshocked inte- 
rior to ~ 10^ cm~^ on the bright ring of emis- 
sion, indicative of rapid cooling of material passing 
through the shock front, or significant pre-shock 
density inhomogeneity. 

3. The ionization state increases rapidly as the ejecta 
encounter the reverse shock, with low ionization 
species dominating inside the shock, moderate ion- 
ization species strongest just outside the shock, and 
very high ionization material seen in X-rays at the 
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TABLE 2 
Neon Crescent Line Fluxes 



Name 


Position 


Area 


^[Nell] 


^[Ne III] 


^[NeV] ^ 


MNe^ 




J2000 


pc^ 




10-15 W/m2 




10-5 Mq 


Ne South 


23:23:31.53 +58:47:22.6 


0.33 


15.9 + 0.04 


5.59 + 0.06 


<0.09^ 


8.9 + 0.4 


Ne North 


23:23:25.88 +58:50:07.2 


0.21 


16.9 + 0.04 


5.62 + 0.05 


0.735 + 0.02 


9.7 + 0.5 



^ Excluding neutral neon.^ Corrected for Ay = 5 magnitudes of extinction. Statistical uncertainties only. Additional 
systematic uncertainties apply (~10% for [Neii],[Neiii]; ~25% for [Nev]).^ 3-cr upper limit 



greatest distances. These variations in post-shock 
history may depend on the density structure of the 
inhomogenous ejecta. 

4. Two symmetrically arranged, compact, crescent- 
shaped clumps of material with highly enhanced 
neon abundance (> 10 x the remnant average) lie 
along the projected direction of the "kick vector" 
separating the compact object from the expansion 
center of the remnant, tracing a new axis of ex- 
plosion asymmetry in the ejecta. Since neon is ex- 
pected to arrive at the reverse shock before heav- 
ier elements like Si and S, these could represent 
directions of reduced initial velocity, intermediate 
between the bright rim, and the (still) quiescent 
material flowing towards the front and back faces 
of the remnant, which has yet to reach the reverse 
shock. 

5. Since the ionization state and density vary steeply 



with distance from the expansion center, and have 
such a strong impact on the observed strength of 
individual emission lines, the apparent macroscopic 
mixing of elements in Cas A may result from vari- 
ations in initial outflow velocity rather than bulk 
turbulent redistribution of the ejecta strata. The 
symmetry in mass and outflow direction of the neon 
crescent regions strongly supports this interpreta- 
tion. 
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